Quantum-transport theory for semiconductor nanostructures: 
A density-matrix formulation 



in 
o 
o 

(N 
Or 

c/3 

(N 



O 

CO 



Rita C. Iotti, Emanuele Ciancio, and Fausto RossQ 
Istituto Nazionale per la Fisica della Materia (INFM) and Dipartimento di Fisica, 
Politecnico di Torino, Corso Duca degli Abruzzi 24, 10129 Torino, Italy 
(Dated: February 2, 2008) 

A general density-matrix formulation of quantum-transport phenomena in semiconductor nanos- 
tructures is presented. More specifically, contrary to the conventional single-particle correlation 
expansion, we shall investigate separately the effects of the adiabatic or Markov limit and of the 
reduction procedure. Our fully operatorial approach allows us to better identify the general prop- 
erties of the scattering superoperators entering our effective quantum-transport theory at various 
description levels, e.g., N electrons-plus-quasiparticles, N electrons only, and single-particle picture. 
In addition to coherent transport phenomena characterizing the transient response of the system, 
the proposed theoretical description allows to study scattering induced phase coherence in steady- 
state conditions. As prototypical example, we shall investigate polaronic effects in strongly biased 
semiconductor superlattices. 

PACS numbers: 72.10.-d, 85.35.-p, 73.63.-b 
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I. INTRODUCTION 

Recent developments in nanoscience/technology 
pushes device miniaturization toward limits where the 
traditional semiclassical transport treatments^ can 
no longer be employed, and more rigorous quantum- 
transport approaches are imperative^ However, in 
spite of the quantum-mechanical nature of carrier 
dynamics in the core region of typical nanostructured 
devices — like semiconductor superlattices, double- 
barrier structures, and quantum dots — the overall 
behavior of such quantum systems is often the result 
of a non-trivial interplay between phase coherence and 
energy relaxation/dephasing£ It follows that a proper 
treatment of such novel nanoscale devices requires 
a theoretical modelling able to properly account for 
both coherent and incoherent — i.e., phase-breaking — 
processes on the same footing within a many-body 
picture. 

More precisely, the idealized behavior of a so-called 
"quantum device"— is usually described via the elemen- 
tary physical picture of the square- well potential and/or 
in terms of a simple quantum-mechanical n-level sys- 
tem. For a quantitative investigation of state-of-the-art 
quantum optoelectronic devices, however, two features 
strongly influence and modify such simplified scenario: 
(i) the intrinsic many-body nature of the carrier sys- 
tem under investigation, and (ii) the potential coupling 
of the electronic subsystem of interest with a variety of 
interaction mechanisms, including the presence of spatial 
boundaries ££il 

The wide family of so-called quantum devices can be 
divided into two main classes: a first one grouping semi- 
conductor devices characterized by a genuine quantum- 
mechanical behavior of their carrier subsystem, and a 
second one which comprises low-dimensional nanostruc- 
tures whose transport dynamics may be safely treated 
within the semiclassical picture. 



Devices within the first class — characterized by a 
weak coupling of the carrier subsystem with the host 
material — are natural candidates for the implementa- 
tion of quantum information/computation processing. 8 
These include, in particular, semiconductor quantum- 
dot structures^ for which all-optical implementations 
have been recently proposed. 10 : 11 : 12 In this case, the pure 
quantum-mechanical carrier dynamics is only weakly dis- 
turbed by decoherence processes; therefore, the latter are 
usually described in terms of extremely simplified mod- 
els. 

Conversely, quantum devices in the second class — in 
spite of their partially discrete energy spectrum due to 
spatial carrier confinement — exhibit a carrier dynamics 
which can be still described via a semiclassical scatter- 
ing picture. Such optoelectronic nanostructured devices 
include multi-quantum- well and superlattice structures, 
like quantum-cascade lasers (QCLs)iiiii These systems 
are characterized by a strong interplay between coherent 
dynamics and energy-relaxation/dephasing processes; it 
follows that for a quantitative description of such non- 
trivial coherence/dissipation coupling the latter need to 
be treated via fully microscopic models. 

In this paper we shall primarily focus on this sec- 
ond class of quantum devices, providing a comprehen- 
sive microscopic theory of charge transport in semicon- 
ductor nanostructures based on the well-known density- 
matrix approach. It is worth mentioning that an alter- 
native approach, equivalent to the density-matrix formal- 
ism employed in this paper, is given by the nonequilib- 
rium Green's function technique; the latter can be re- 
garded as an extension of the well-known equilibrium or 
zero-temperature Green's function theory to nonequilib- 
rium regimes, introduced in the 1960s by Kadanoff and 
Baymi^ and Keldysh>i& An introduction to the theory 
of nonequilibrium Green's functions with applications to 
many problems in transport and optics of semiconduc- 
tors can be found in the book by Haug and Jauho^i 
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By employing — and further developing and extending — 
such nonequilibrium Green's function formalism, a num- 
ber of groups have proposed efficient quantum-transport 
treatments for the study of various semiconductor nanos- 
tructuresi^ as well as of modern micro/optoelectronic de- 
vices. 

Within the general density-matrix formalism two 
different strategies are commonly employed: (i) the 
quantum- kinetic treatment^ and (ii) the description 
based on the Liouville-von Neumann equation^ 

The primary goal of a quantum-kinetic theory is to 
evaluate the temporal evolution of a reduced set of single- 
or few-particle quantities directly related to the electro- 
optical phenomenon under investigation, the so-called ki- 
netic variables of the system. However, due to the many- 
body nature of the problem, an exact solution in general 
is not possible; it follows that for a detailed understand- 
ing realistic semiconductor models have to be considered, 
which then can only be treated approximately. Within 
the kinetic-theory approach one starts directly with the 
equations of motion for the single-particle density ma- 
trix. Due to the many-body nature of the problem, the 
resulting set of equations of motion is not closed; instead, 
it constitutes the starting point of an infinite hierarchy of 
higher-order density matrices. Besides differences related 
to the quantum statistics of the quasiparticles involved, 
this is equivalent to the BBGKY hierarchy in classical gas 
dynamics^ The central approximation in this formalism 
is the truncation of the hierarchy. This can be based on 
different physical pictures. A common approach is to use 
the argument that correlations involving an increasing 
number of particles will become less and less important. 3 
An alternative quantum-kinetic scheme — based on an 
expansion in powers of the exciting laser field — has been 
introduced by Axt and Stahl, the so-called "dynamics 
controlled truncation" (DCT)^ 

Within the treatment based on the Liouville-von Neu- 
mann equation, the starting point is the equation of mo- 
tion for the global density-matrix operator, describing 
many electron plus various quasiparticle excitations. The 
physical quantities of interest for the electronic subsys- 
tem are then typically derived via a suitable "reduction 
procedure" , aimed at tracing out non-relevant degrees of 
freedom. Contrary to the kinetic theory, this approach 
has allowed for a fully quantum-mechanical treatment 
of high-field transport in semiconductors, 23 thus over- 
coming some of the basic limitations of conventional ki- 
netic treatments, e.g., the completed-collision limit and 
the Markov approximation. 

Primary goal of the present paper is to discuss in very 
general terms the physical properties and validity limits 
of the so-called "adiabatic" or Markov approximation. 
Within the traditional semiclassical or Boltzmann the- 
ory, this approximation is typically introduced together 
with the so-called diagonal approximation, i.e., the ne- 
glect of non-diagonal density-matrix elements. However, 
as described in Ref. 1241 the Markov limit can be also per- 
formed within a fully non-diagonal density-matrix treat- 



ment of the problem; this leads to the introduction of 
generalized in- and out-scattering superoperators, whose 
general properties and physical interpretation are not 
straightforward. In particular, it is imperative to under- 
stand if — and under which conditions — the adiabatic 
or Markov approximation preserves the positive-definite 
character of our reduced density matrix; indeed, this dis- 
tinguished property is generally lost within the quantum- 
kinetic approaches previously mentionedA^i To this end, 
starting from the Liouville-von Neumann-equation ap- 
proach, we shall propose a very general treatment of the 
Markov approximation. More specifically, contrary to 
the conventional single-particle correlation expansion of 
the kinetic theory, we shall investigate separately the ef- 
fects of the Markov limit and of the reduction procedure. 
Our fully operatorial approach will allow us to better 
identify the general properties of the scattering super- 
operators entering our effective quantum-transport the- 
ory at various description levels, e.g., N electrons-plus- 
quasiparticles, N electrons only, and single-particle pic- 
ture. In addition to coherent transport phenomena char- 
acterizing the transient response of the system, the pro- 
posed theoretical description allows to study scattering 
induced phase coherence in steady-state conditions. In 
particular, based on the proposed approach we shall con- 
sider — as prototypical example — polaronic coherence in 
strongly biased semiconductor superlattices. 

The paper is organized as follows. In Sect. [H] we 
present the proposed theoretical approach: after speci- 
fying the physical system under investigation fSect. lHA'l) 
and recalling the fundamentals of the density-matrix for- 
malism (Sect. Ill B|) . we shall introduce the Markov ap- 
proximation (Sect. Ill C|) and derive the explicit form of 
the scattering superoperators (Sect. lHD|) as well as their 
semiclassical counterparts fSect. UlE}) : we shall then dis- 
cuss the so-called reduction procedure fSect. lHF|l and the 
single-particle description (Sect. Ill G|) for carrier-carrier 
as well as carrier-quasiparticle interactions. In Sect. IIIII 
we shall address the general problem of scattering- 
induced phase coherence; in particular, we shall present 
a few simulated experiments concerning polaronic coher- 
ence in semiconductor superlattices. Finally, in Sect. |V] 
we shall summarize and draw some conclusions. 



II. THEORETICAL BACKGROUND 

A. Physical system 

To provide a general formulation of quantum charge 
transport in semiconductor nanostructures, let us con- 
sider a generic carrier gas within a semiconductor crys- 
tal in the presence of electromagnetic fields. The corre- 
sponding Hamiltonian can be schematically written as 

H = H a +H' . (1) 
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The first term, 



£ q fe q fe q 



(2) 



is the sum of the free-carrier and free-quasiparticle 
Hamiltonians, where the Fermionic operators & a (c a ) 
denote creation (destruction) of a carrier in the single- 
particle state a (with energy e Q ), while the Bosonic op- 
erators & q (6 q ) denote creation (destruction) of a generic 
quasiparticle excitation with wavevector q and energy e q , 
i.e., phonons, photons, plasmons, etc. 

The second term, H', is the sum of all possible 
carrier-carrier as well as carrier-quasiparticle interac- 
tion Hamiltonians, i.e., carrier-phonon, carrier-photon, 
carrier-plasmon, etc. 

The noninteracting carrier-plus-quasiparticle basis sta- 
tes are given by the eigenstates of H : the generic eigen- 
state | A) = |{n q }) is the tensor product of non- 

interacting carrier and quasiparticle states correspond- 
ing, respectively, to the occupation numbers {n a } and 
{n q }, while the noninteracting energy spectrum e\ = 
J2 a e a n a + S q e q n q is the sum of the total carrier and 
quasiparticle energies. 

The interaction Hamiltonian H' cannot in general be 
treated exactly. A typical approach consists of regarding 
it as a perturbation acting on the noninteracting carrier- 
plus-quasiparticle states {|A)}. In this context, the basic 
ingredients are the matrix elements of H ' within our non- 
interacting basis states: H' xx , = (X\H'\X'}. 



B. Density-matrix formalism 

In view of the huge number of degrees of freedom 
({a}, {q}) involved in the microscopic treatment of any 
solid-state system, a statistical description of the prob- 
lem is imperative. As we shall see (in Sect. Ill F|) . this will 
result in a suitable statistical average over "non-relevant" 
degrees of freedom. 

Given a physical quantity A — described by the op- 
erator A — its quantum plus statistical average value is 
given by 
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^\A\^)=tx{Ap) , 



where 



p=\il>M 



(3) 



(4) 



is the so-called density-matrix operator. The latter is de- 
fined as statistical average of the projection operator cor- 
responding to the generic state vector of the system, 
and can then be regarded as the statistical generalization 
of the quantum-mechanical concept of state vector. 

Starting from the global Schrodinger equation describ- 
ing our interacting carrier-plus-quasiparticle many-body 
system, the following Liouville-von Neumann equation 



of motion for the density-matrix operator can be readily 
obtained: 



— = L (p) = — 
dt ih 



H,p 



(5) 



where C is usually referred to as Liouville superoperator. 
Equation © can be regarded as the statistical general- 
ization of the Schrodinger equation; Its exact solution is 
given by 

p(t) = e c ^p(t Q ) = U(t - t )p(t )U\t - t ) , (6) 
where 



U(t-t ) = e 



H(t-t a ) 



(7) 



is the evolution operator corresponding to the total 
Hamiltonian H in Such exact solution corre- 

sponds to a fully quantum-mechanical unitary evolution 
of the whole many-body system, i.e., no energy relax- 
ation/dephasing. Indeed, it is easy to verify that the 
total quantum entropy 



S = -fc B tr{plog/5} 



(8) 



is not affected by the unitary transformation in Q. 25 
As anticipated, the total many-body Hamiltonian in 
cannot be treated exactly. Aim of a quantum- 
transport theory is to derive effective equations describ- 
ing the carrier subsystem of interest within some approx- 
imation scheme; this is typically realized via the follow- 
ing two basic steps: first an adiabatic decoupling between 
the different time-scales induced by H and H' — called 
Markov limit — and then a projection of the global sys- 
tem dynamics over a subsystem of interest via the intro- 
duction of a so-called reduced density-matrix operator. 



C. The adiabatic or Markov approximation 

Starting from the separation H = H a + H' in (JIJ, the 
Liouville-von Neumann equation (J5J can be written as 



dp dp 
~dl ~ ~dt 



H, 



dp 
~dl 



H' 



(9) 



where the two contributions describe, respectively, the 
time evolution induced by the noninteracting Hamilto- 
nian H and by the interaction term H' . 

The first contribution can be treated exactly within the 
standard interaction scheme. Indeed, it is easy to show 
that the time evolution of the density-matrix operator in 
the interaction picture, 



p l = Ul{t-t )pU o (t-t ) 



is simply given by 



dp 1 
~dt 



(10) 



(11) 
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where U (t — to) is the evolution operator corresponding 
to the noninteracting Hamiltonian H , and 

H i {t) = U${t-t ) I ^-U o {t-t Q ) (12) 

denotes the Hamiltonian H' in units of h within the in- 
teraction picture. 

The key idea beyond any perturbation approach is that 
the effect of the interaction Hamiltonian W is "small" 
compared to the free evolution dictated by the noninter- 
acting Hamiltonian H . More precisely, the interaction 



J 



matrix elements H' xx , are smaller than the typical energy 
difference e\ — ey . 

Following this spirit, by formally integrating Eq. I|ll|) 
from to to the current time i, we get: 



p\t)=f?{t )-i / dt' 
Jto 



(13) 



By inserting the above formal solution for p l (t) on 
the right-hand side of Eq. I|ll|) we obtain an integro- 
differentia! equation of the form: 



dt' 



r 



(14) 



We stress that so far no approximation has been intro- 
duced: Equations ifTTJl. Q13|). and (|H|l are all fully equiv- 
alent, we have just isolated the first-order contribution 
from the full time evolution in Eq. (|llfl . It is then clear 



that, by iteratively substituting Eq. (|13|) into itself, the 
above procedure can be extended to any perturbation 
order. This leads to the well-known Ncuman series: 



! (<) = p ! (M + ]LHr / dt * / dt * 



to Jto 



dt n [H\t l ),[HXt 2 ),...[H\t n ),p\t a )} 



(15) 



r 



The latter constitutes the starting point of the quantum 
Monte Carlo method for the study of charge-transport 
phenomena in semiconductors^ 

In order to introduce the so-called adiabatic or Markov 
approximation, let us now focus on the time integral in 
Eq. lfl4T) . Here, the two quantities to be integrated over 
t' are the interaction Hamiltonian TV and the density- 
matrix operator p l . In the spirit of the perturbation ap- 
proach previously recalled, the time variation of p 1 can 
be considered adiabatically slow compared to that of the 
Hamiltonian TC 1 within the interaction picture; indeed, 
the latter will exhibit rapid oscillations due to the non- 
interacting unitary transformation U a . As a result, the 
density-matrix operator p l can be taken out of the time 
integral and evaluated at the current time t. 

Within such adiabatic limit we get the following effec- 
tive Liouville-von Neumann equation: 



dt' 
with 



p\t) = -i n l (t),p l (t ) - n\t), ic\t),p l (t) 



JC{t) = / H l {t')dt' 

Jt 



(16) 
(17) 



The above equation has still the double-commutator 
structure in l|14|) but it is now local in time. 

Going back to the original Schrodinger picture, we fi- 
nally get: 



dp 
~dl 



C(p) + C 



where 



with 



1 r 



H a ,P 



ic, p 



1 r 

if i 



H Q ,p 



K. 



dt'U {t-t')HUt{t-t') , 



(18) 

r(p) , 

(19) 
(20) 



is the effective Liouville superoperator within our approx- 
imation scheme, and 



C(t) 



H,U o (t-t Q )p(t Q )Ul(t-t ) 



(21) 



The time-dependent operator C in (|21|l describes how 
the quantum-correlation effects at the initial time to 
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propagate to the current time t; indeed, combining 
Eqs. © and H21|) . the latter can be rewritten as: 



C(i) 



n,s(t-t )p(t)s^{t-t ) 



-i\nrp % {t)\ , 

(22) 

where S(t — t ) = U (t — t a )U\t — t ) is the uni- 
tary transformation connecting the time evolution of 
the density-matrix operator in the Schrodinger and in- 
teraction pictures. This clearly shows that the ini- 



J 



tial quantum-mechanical correlations propagate from to 
to t via the interaction-free dynamics described by the 
density-matrix operator written in the interaction pic- 
ture. As we shall see, the above quantum-correlation 
operator is responsible for a number of purely quantum- 
mechanical phenomena, like Hartree-Fock single-particle 
renormalizations and coherent phonon effects^ 

The general solution of Eq. I|18f) is of the form: 



ffl = T 



3 tn 



dt' 



P C(t")dt" 



C{t')dt' , 



(23) 



where T [. . .] is the usual time- or chronological-ordering 
operatorii 5 ^ 

At this point a few comments are in order. So far, the 
only approximation introduced in our theoretical descrip- 
tion is the adiabatic decoupling between free carrier evo- 
lution and various many-body interactions; this leads to 
a significant modification of the system dynamics: while 
the exact quantum-mechanical evolution in corre- 
sponds to a fully reversible and isoentropic unitary trans- 
formation, the instantaneous double-commutator struc- 
ture in (|19[) describes, in general, a non-reversible (i.e., 
non unitary) dynamics [see Eq. (|23[) ] characterized by 
energy relaxation and dephasing; it follows that the sys- 
tem quantum entropy in JSJ is no more a constant. At 
this level of description this behavior is totally nonphys- 
ical, clearly showing the potential failure and intrinsic 
limitations of the Markov approximation. However, as 
discussed below (see Sect. Ill F|l . the Markov limit pre- 
viously introduced is usually employed together with a 
reduced description of the system, for which such irre- 
versible dynamics is physically justified. 

Let us finally focus on the nature of the effective Li- 
ouville superoperator in (|19|) . As stressed before, this 
is the sum of a single-commutator term plus a double- 
commutator contribution. In the absence of carrier- 
carrier as well as carrier-quasiparticle interactions, i.e., 
H' = 0, the second term vanishes and the system un- 
dergoes a reversible unitary transformation induced by 
the single-commutator term, which preserves the trace 
and the positive character of our density-matrix oper- 
ator p. In contrast, the perturbation Hamiltonian H' 
within the Markov limit previously introduced will in- 
duce, in general, a non-unitary evolution. Since any ef- 
fective Liouville superoperator should describe correctly 
the time evolution of p and since the latter, by defini- 
tion, needs to be trace-invariant and positive-definite at 
any time, it is important to determine if — and under 
which conditions — the superoperator C fulfills this two 
basic requirements. 

As far as the first issue is concerned, recalling that the 



trace of a commutator is always equal to zero and taking 
the trace of Eq. (|18fl . it is easy to verify that the time- 
derivative of the trace of p is equal to zero, i.e., that our 
effective dynamics is trace-preserving. 

Let us now discuss the possible positive-definite char- 
acter of p. In general, our effective Liouville superoper- 
ator does not ensure that for any initial condition the 
density-matrix operator will be positive-definite at any 
time. Indeed, it is possible to show that the double- 
commutator structure in l|19l) can be rewritten in terms 
of a single-commutator structure (renormalizing the free 
Hamiltonian H ) and of double commutators of the form: 



L(p) = - A, A,p 



(24) 



Each of the latter represents a particular case of the 
so called Lindblad superoperators, which are known to 
describe completely-positive (CP) maps, thus preserv- 
ing the positive character of our density-matrix operator. 
However, our Liouville superoperator can be written in 
terms of the difference of the Lindblad superoperators in 
(|24|l . which in general is not Lindblad-like. 

Since our primary goal is the investigation of quantum- 
transport phenomena, we shall focus on the steady-state 
solution of Eq. I|18|) ■ It is easy to verify that the identity 
operator, properly normalized, p(t) oc T is the station- 
ary solution we are looking for. A closer inspection of 
Eq. lll'/il) reveals that for any positive-definite and un- 
corrected initial state (C — 0) and for a Liouville su- 
peroperator C with a non-positive eigenvalue spectrum, 
in the limit t — > oo the density matrix p will reduce to 
the identity operator previously mentioned. As antici- 
pated, this clearly shows that within such approximation 
scheme our effective dynamics describes a sort of decoher- 
ence/dephasing, since possible non-diagonal terms of the 
density matrix will vanish on the long-time scale. This 
is again an artefact of the Markov limit. 

Let us finally discuss the physical meaning of the 
steady-state solution p(t — > oo) oc 2. Within our non- 
interacting carrier-plus-quasiparticle basis A we have: 
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PAiA 2 (i — * °°) ^AiA 2 - This tells us that, physically 
speaking, the steady-state solution of our transport equa- 
tion corresponds to an equally-probable population of all 
the microscopic states A without any interstate quan- 
tum coherence (/CAi^a 2 = 0). This scenario is typical of 
the present global (carrier + quasiparticle) description; 
in contrast, within a reduced description (of the carrier 
subsystem only) the steady-state solution differs from the 
identity operator, since in this case the trace over non- 
relevant degrees of freedom will translate into a thermal 



weight over our electronic states (see Sect. Ill Ft 



D. Generalized scattering superoperator 

Let us now evaluate the explicit form of the scatter- 
ing superoperator T in i|2(J|) . The effective Liouville-von 
Neumann equation in (|18|l can be easily rewritten in the 
noninteracting-states basis {|A)} previously introduced: 



J 



dt 



CAi — £A £ 

if) 



px t x 2 + C\ t \ 2 + Fai x 2 ,x' l x' 2 Px' 1 x 2 



(25) 



In order to derive the explicit form of the superopera- commutator in H2Q(I : 



tor matrix elements T 



Ai A 2 , A^ X' 2 j 



let us expand the double 



FAi X2XK K — K^' \fciP , , — (^AiA'PA'jA'^A' A 2 + ^AiA'.PA'.A' Wa'A 2 

L L J J Ai A 2 

y ± x 2 ^'1^2 

- X] (Wa 1 A^A' 1 A^PA 2 A 2 + PAiA'^A'jA^A^Aa) ■ 



(26) 



As we can see, the scattering operator T can be written as 
the difference of the following "in-" and "out-scattering" 
terms: 

^Ai A 2 ,A' 1 Aj = ^A 1 A' 1 /C A ^A 2 + ICXiX'jHx'^ (27) 



r AiA 2 ,AiA^ - ^2'Hx 1 x>>ICx"x> 1 h 2 x 2 + ^2h 1 x> 1 fcx> 2 x>>'H 



A"A 2 



(28) 



r 



For the case of a time-independent perturbation H' = with 
KH, the operator fC can be rewritten as: 



£= ! ° drU {r)HUt{T) 
Jo 



(29) 



e dr = VX, 



A'A 



(32) 



Taking into account that within the A-representation the 
noninteracting evolution operator U is simply given by 

U xx '(T) = e^S xx , (30) 
the matrix elements of the operator K. in l(2*9"]l result to By inserting the above result into Eqs. lfT7j) - lf28|) and re- 



be: 



ICxx' — ^'Hxx'T^xx 1 



(31) 



calling that 7i 



XX' 



"HX/x, we finally obtain: 
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Ai A2 ,A1 XL 



2tt / 



AiA, A9A: 



Aq A' 



(33) 



pout 

Ai A2 jAj^ A 2 



2tt 



E 

A" 



ttI* 



(34) 



In general, the scattering superoperator T is a function 
of time; however, in the limit t Q — > — 00 — also called 
"completed-collision limit"^^ — the function V in (1321) 
becomes time-independent: 

^aT = e^^dr . (35) 

It follows that in this limit the operator JC as well as the 
superoperators T and C are also time-independent. In 
this case there is no need for the time-ordering operator 
T in 123(1 . Moreover, the real part of the function T>~°° in 
((35(1 gives the well-known energy-conserving Dirac delta 
function, i.e., 

= ^5(e x -e x ,)+K X y , (36) 

while the imaginary part — denoted by IZw — describes, 
in general, energy-renormalization effects. Within the 
validity limits of the present Markov treatment, such 
renormalization effects can be safely neglected: if, as re- 
quested, the perturbation Hamiltonian is small compared 
to the noninteracting one, then the resulting energy-level 
renormalization is small compared to the noninteracting 
energy levels e\. 

At this point a few comments on the evaluation of 
the time integral in 1(35(1 are in order. Indeed, it is well 
known that the limit to — > — 00 needs to be performed 
properly; more specifically, this is realized by adding to 
the energy difference an infinitesimally small imaginary 
part, which ensures the convergence of the time inte- 
gration. A qualitative — but not rigorous — interpreta- 
tion of such mathematical prescription is based on the 
so-called "adiabatic switching-on" procedures^ The idea 
is that, starting from t = — 00, the interaction mecha- 
nism/Hamiltonian is slowly or adiabatically switched on. 
By employing the noncquilibrium Green's function for- 
malism, it is possible to showiZ. that such imaginary part 
is not an artificial ingredient: it corresponds to the imag- 
inary part of the electron self-energy, thus describing the 
finite life-time of our electronic states due to all relevant 
interaction mechanisms. A proper account of such ef- 
fect — not relevant in the present discussion — leads to 
apparent violations of the energy-conserving transitions 
predicted by the Dirac delta function in ((36(1 . the so- 
called "collisional broadening" . 

As previously recalled, it is imperative to establish if 
— and under which conditions — the scattering superop- 



erator in ((26(1 preserves the positive-definite nature of our 
density-matrix operator p. As discussed in App. 

(i) contrary to the semiclassical or Boltzmann dynam- 
ics (see Sect. Ill E|) . the effective Liouville superop- 
erator previously identified does not correspond to 
a so-called "CP map", i.e., it does not preserve, 
in general, the positive-definite character of our 
density-matrix operator; 

(ii) its eigenvalue spectrum, i.e., 

£ (p) = Ap , (37) 

always contains the A = eigenvalue, which corre- 
sponds to the steady-state transport solution; 

(iii) in the "small-coupling limit" it is possible to show 
that the steady-state density-matrix operator - 
corresponding to the A = eigenvalue — is always 
positive definite, i.e., 

P A =° = 1>/3IWI > ( 38 ) 
A 

where the basis states |/3) are the eigenvectors of 
p A=0 , and Pa is a (non- negative) probability distri- 
bution. 



E. Semiclassical limit 

The well-known semiclassical or Boltzmann trans- 
port theorjii can be easily derived from the quantum- 
transport formulation presented so far, by introducing 
the so-called diagonal or semiclassical limit. The latter 
corresponds to neglecting all non-diagonal density-matrix 
elements (and therefore any quantum-mechanical phase 
coherence between the generic states Ai and A2), i.e., 

P\ x \2 = /a!<5aiA 2 , (39) 

where the diagonal elements f\ describe the semiclassical 
distribution function over the noninteracting basis states 
A. 

By introducing the above semiclassical density matrix 
into Eq. ((25|) for the diagonal elements (A = Ai = A2), 
and inserting the explicit form of the elements T\\ of the 
scattering operator in the limit to — > —00 [see Eq. (|35(lj. 



8 



we finally obtain the usual form of the Boltzmann trans- 
port equation written in our basis states: 



^ = £(Paa'/a'-JWa) 

A' 



where 



P: 



AA' 



Py 



A' A 



2tt 

T 



H' x , x \ 2 5(e x , 



ex) 



(40) 



(41) 



are the semiclassical scattering rates given by the well- 
known Fermi's golden rule. In addition to the square 
of the interaction matrix element Hy X , they contain the 
energy-conserving Dirac delta function: 



S {tx> - e\) 



V 



AA' 



^A'A 



1 

2^ 



Ae x ,-ex)T dT 



(42) 

Within the semiclassical limit the free-rotation term in 
(|25|l vanishes, and the same applies to the quantum- 
correlation contributions C\\. 

Our analysis shows that the quantum-transport equa- 
tion in l)25|) can be regarded as the quantum-mechanical 
generalization of the Boltzmann equation in (|40|) . In- 
deed, the in- and out-scattering superoperators in Ij26(l 
are the quantum-mechanical generalizations of the stan- 
dard in- and out-scattering terms entering the Boltzmann 
collision operator in (|40|) . 

As a confirmation of the fact that the Markov approxi- 
mation leads to a totally nonphysical non-reversible (i.e., 
non-unitary) system evolution, it is possible to show that 
the system entropy S in Eq. ((SJ) is a non-decreasing func- 
tion of time. 

We stress that, contrary to the usual semiclassical 
transport theory, the Boltzmann-like equation in 14011 
describes a scattering dynamics within the whole A = 
{n Q },{n q } space, i.e., the generic scattering probabil- 
ity P\y in 141|) describes a transition from the state 
{n a }, {riq} to the state {n a >}, {n q <}. In other words, 
so far no reduction procedure to the a-subsystem has 
been performed; it follows that for a given transition 
of the a subsystem ({n a } — > {n a >}), a corresponding 
transition ({n q } — > {n q /}) of the quasiparticle subsys- 
tem will also take place. This explains why, contrary to 
the usual Boltzmann theory, the scattering probabilities 
in 1411) are symmetric, i.e., invariant under time rever- 
sal: P\\r = Px'X'i moreover, the Dirac delta function in 
(|41fl leads to the conservation of the total energy of the 
system. 

A second important remark is that, contrary to the 
non-diagonal density-matrix description previously intro- 
duced, the Markov limit combined with the semiclassical 
or diagonal approximation in l|39l) ensures that at any 
time t our semiclassical distribution function f\ is always 
positive-definite (see App. IA"|). 

Let us finally discuss the steady-state solution of 
the Boltzmann transport equation in H40|) . From the 
detailed-balance principle, i.e., 



and considering that the semiclassical scattering rates in 
(|4*T1) are symmetric (Pxx< = Px'x), we get: 



fx 



Pxx> 

p. 



= 1 



A' A 



(44) 



Exactly as in the quantum-mechanical case, the steady- 
state solution corresponds to a uniform distribution over 
the noninteracting carrier-plus-quasiparticle states: fx oc 
S\\ . As discussed in Sect. Ill Fl this is not the case when 
our kinetic description is reduced to the a subsystem 
only. 



F. Reduced description 

As discussed in Sect. Ill Bl the average value of any 
given physical quantity A can be easily expressed in terms 
of the density-matrix operator p according to Eq. yj. 
In the study of charge-transport phenomena in semicon- 
ductor nanostructures, most of the physical quantities 
of interest depend on the electronic-subsystem coordi- 
nates a only (carrier drift velocity, total electronic en- 
ergy, carrier-carrier correlation function, etc.), i.e., 

A M, {»,);{%:), {r V } = A {n a }.{n al } 5 {n^},{n w } ■ (45) 

In this case it is convenient to write 
A = ^2 A XX'PX'X = ^2 ^Kl.lv }P{n a , },{«<*} ' 

AA' {n Q } , { n Q / } 

(46) 

where 

p {n a },{n a >} = H ^W. {«<,};{»='}.{«<,} ( 4? ) 
{«„} 

is the so-called reduced or electronic density matrix. 
Equation (|47|l can be also written in an operatorial form 



p c = ti{p} 



{™q} 



(48) 



which shows that the electronic density-matrix operator 
p c is obtained by performing a trace operation over the 
quasi-particle coordinates q. Since p c is the only quan- 
tity entering the evaluation of the average value in I4(ill , 
it is desirable to derive a corresponding equation of mo- 
tion for the reduced density-matrix operator. Combining 
Eqs. (HHJ and gHl we get: 



f=tr{£(p)| +tr{(7} 

dt I J {n q } I J {n q } 



(49) 



Pxx'fx' = Pyxfx 



(43) 



In general, the trace over the quasiparticle coordinates 
does not commute with the Liouville superoperator L in 
(|19fl . which does not allow to obtain a closed equation 
of motion for the reduced density-matrix operator p c . 
This clearly does not apply when the interaction Hamil- 
tonian is a function of the carrier coordinates only, e.g., 
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carrier-carrier, carrier-impurity, etc. In contrast, in the 
presence of carrier-quasiparticle coupling additional ap- 
proximations are needed. In order to get a closed equa- 
tion of motion for the reduced density-matrix operator, 
the typical assumption is to consider the quasiparticle 
subsystem as characterized by a huge number of degrees 
of freedom (compared to the subsystem a). In other 
words this amounts to say that the q subsystem has an 
infinitely high heat capacity, i.e., it behaves as a thermal 
bath; this allows to consider the quasiparticle subsystem 
always in thermal equilibrium, i.e., not significantly per- 
turbed by the carrier subsystem a. Within such approx- 
imation scheme, the global (a + q) density-matrix oper- 
ator p can be written as the product of the equilibrium 
density-matrix operator for the quasiparticle subsystem 
p qp and the reduced density-matrix operator p c : 

p = p c ®p^, p™ = r ^ 1 ■ (50) 

tr <^ e "b t I 



The corresponding matrix elements within our basis 
states A = {n a }, {n q } are then given by: 

PXX' = P{n a }{n a ,}f 9 { n,} S {n^{n^} (51) 

with 

_!W> 

m } - - — w ■ (^) 

£{n,} e *» T 

By inserting these density-matrix elements into Eq. Ij25(l 
and performing the trace over the quasiparticle coordi- 
nates, it is easy to get the following effective Liouville-von 
Neumann equation for the reduced density matrix p c : 



{n Q1 }{"c, 2 } 

dt 



i% ' , NHsl + G K}W + 2^ 1 KlK ! } ! {vHv}''{v}{vl 

{»«'}{•»„'} 



(53) 



with 



{n ai }{n a2 },{n, }{n, } 



^ -^{ n £»l}{™q}'{™a 2 }{ Il q};{™Q ! ' 1 }{™q'}d™ Q ^}{» l q '}^{ri cl /} 
{™q}{"q'} 



(54) 



and 



tonian H' can be written as: 



C {n ai }{n a2 } - J! ^{««lH«<J>{»<* 2 }{»<>} • 

{»q} 



(55) 



By denoting with r c the effective scattering superopera- 
tor defined by the matrix elements in (|54|l — acting on the 
a Hilbert subspace only — the new effective Liouville su- 
peroperator (i.e., traced over the q coordinates) is given 
by: 



t {p c 



I 

if l 



m 



p 



+ r c (p c ) 



(56) 



The latter, however, does not contain the double- 
commutator structure previously discussed; This aspect 
will be more extensively addressed in the following. 

For all relevant carrier-quasiparticle interaction mech- 
anisms in semiconductor nanostructures — e.g., carrier- 
phonon, carrier-plasmon, etc. — the perturbation Hamil- 



h' = nn = n 



(57) 



Here 7i q = 7iL q are electronic operators (parameterized 
by the quasiparticle wavevector q) acting on the a sub- 
system only. The two terms in (|57|) — corresponding to 
quasiparticle destruction and creation — describe quasi- 
particle absorption and emission processes. 

Let us consider again the definition of the effective 
carrier-quasiparticle scattering superoperator T c in l|54|) 
written in operatorial form, i.e., 



r c (p c ) = -tv{\n, \ic,p c H]} , 

U L J J > {«q} 



(58) 



where the double-commutator form in l|20l) has been in- 
troduced. 

By inserting into Eqs. I|29|l and (|58|l the explicit form 
of the carrier-quasiparticle Hamiltonian in l|57|l and using 
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= <W> we ob- 
tain an explicit form of the effective carrier-quasiparticle 



the bosonic commutation relations [6 q , 6 q ,j - < - <iq . 



scattering superoperator in ffify . More specifically we 
get: 



r c (P c ) 



with 



E 

q± 









4) 







(59) 



denotes the equilibrium average occupation number for 
the quasiparticle state q. As we can see, for each quasi- 
particle state q we have two contributions (±) describing 
quasiparticle emission and absorption. 

The effective scattering superoperator T c in (|59|l does 
not exhibit the double-commutator structure typical of 
the global description [see Eq. 1(20)1], Indeed, denoting 
with 



Here 



A q = tr 



1 J e fc s T - 1 



(60) 



%i=\ + . ^2 = \ tf) , (62) 



(61) the superoperator T c in l|59|) can be also expressed as: 
I 



r c (p c ) 



E 

q± 



r 



n, 



(63) 



As anticipated, the scattering superoperator involves 
again a double-commutator term, but we have also a 
commutator-anticommutator contribution. 

To better underline the physical role played by the 
above double-commutator versus commutator-anticom- 
mutator contributions, let us recall a simplified model 
usually invoked to qualitatively describe the quantum- 
mechanical evolution of open systems, i.e., subsystems 
interacting with their environment. Within the Schrodin- 
ger picture, the latter are typically treated by adding to 
the system Hamiltonian H an antihermitian part (imag- 
inary potential) y env describing the system-environment 
coupling: 



ih-\1>) = [H + V cm ) |V) 



(64) 



Starting from the above modified Schrodinger equa- 
tion, it is easy to obtain a corresponding version of the 
Liouville-von Neumann equation in 



dp 
~dt 



1 

if i 



H,p 



(65) 



In addition to the commutator-like dynamics typical of a 
closed system, we deal with an anticommutator term, de- 
scribing dissipation induced by the system-environment 
coupling. Indeed, contrary to the closed dynamics in 
(J5J, the latter leads to a non-reversible dynamics. Such 
a simplified model is known to be highly nonphysical, 
since it does not preserve the trace of the density-matrix 
operator; however, it clearly shows how the commuta- 
tor structure is intimately related to a closed evolution, 
while anticommutator terms always describe dissipation 
processes. 



Contrary to the above simplified model, the scattering 
superoperator in (|59|l — in view of its outer-commutator 
structure — is trace-preserving. However, the presence of 
the commutator-anticommutator contribution is a clear 
fingerprint of carrier-quasiparticle dissipation phenom- 
ena leading to genuine energy- relaxation/dephasing pro- 
cesses. 

As far as the correlation term in (|55|l is concerned, it 
is easy to verify that the latter vanishes for the linear- 
coupling carrier-quasiparticle Hamiltonian in l|57|) . 

In summary, within the approximation scheme con- 
sidered so far we get the following effective equation of 
motion for the reduced density-matrix operator p c : 



dp c 
~dt 



t{p c ) 



(66) 



with C c defined in (|56|) . In the limit t — ► — oo, the 
effective Liouville superoperator in (|56|l becomes time- 
independent, and the general solution of the homoge- 
neous equation in l|66|) is of the form: 



P c (t) 



p c (to) 



(67) 



Again, contrary to the isoentropic and fully-reversible 
unitary evolution in ©, the instantaneous double-com- 
mutator plus commutator-anticommutator structures in 
(|63|l describe a non-reversible (i.e., non unitary) dynam- 
ics characterized by energy relaxation and dephasing in- 
duced by the carrier-quasiparticle coupling in (|57Jl . 

We shall now derive the explicit form of the scattering 
superoperator r c within our noninteracting-carrier basis 
{n a }. To this end let us expand the various terms enter- 
ing Eq. JSHJl: 



11 



E T \n al }{n a2 },{n a , }{n a , }P\n a , }{n a , } 

1 2 



E 



q± 



VjK' }{»<} n {" Q ' }{n a2 } 



n c r qT oyq 

^"aJK'}^ HV } n i« Q ' }K> 2 } 
1 12 2 



(68) 



As we can see, also in this case the scattering superoper- 
ator r c can be regarded as the difference of the following 



in- and out-scattering terms: 



J 



r {" Q1 }{n Q2 }{«„;,} _ E + 2 ± 2 J (^JK^'^t'' 



q± 



(69) 



{n 01 }{n Q2 },{n Q , }{r 



: ■} = E ( ^ + 5 ± 5 ) ( ^Hvl^t'Hn^ }*KHV 



{n a «},q± 



r 



(70) 



As for the case of the global (carrier + quasiparticle) 
description presented in Sect. Ill Ul the matrix elements 
JC} n yr ,1 of the operator JC^ in (|fiU|l may be expressed 

in terms of the matrix elements of the operator 7i q : 

JC?*,, x =2irHT x: iZ > ? T \, , (71) 

with 



Again, in the completed-collision limit (to — > -co), 
the real part of I? provides the energy-conserving Dirac 
delta function, while its imaginary part describes carrier- 
quasiparticle energy-renormalization effects. By insert- 
ing this relation into the above in- and out-scattering 



rates and recalling that 7^* w , 



TL r q i r i and 



-pqi 

U {n a }{n a ,} 



1 

2^ 



t-to (< 



"dr . (72) 



q±* 



,q± 



{«d }{»»', } KiUvl {™c« 2 }{"■„/ } 



and 



pout : ^ 

K){»a 2 },{»„;({"^l fi 2 



{n a «}{n a / }^{n a «}{n a / } K 2 H»<,' } 



q± 



(73) 



K»»},q± 



^{loiH™a' }- C> {n //}{n„/ }^{n CI «}{n„/ }^{k.,,}{n ,} 
1 9 r> 



(74) 



where U? 1 = KH^ n j{„ ,\- We stress that the i.e., p- independent. As we shall see, this feature — typical 

above in- and out-scattering superoperators are linear, 
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of the present many-electron description {n a } — will be 
lost in the single-particle picture discussed below (see 
Sect. ESI. 

Let us finally focus on the steady-state solution of the 
quantum-transport equation in l|66() . Contrary to the 
global (carrier + quasiparticle) equation in IjlSjl , the iden- 
tity operator X is no more a solution. Indeed, the latter 
fulfills the double commutator but not the commutator- 
anticommutator structure in l|63|) . Moreover, as we shall 
see in Sect. IIIII the steady-state solution of our effec- 
tive transport equation is in general non-diagonal. We 
also stress that, as for the global (carrier + quasiparti- 
cle) description, in the small-coupling limit it is possible 
to show that the steady-state reduced density matrix is 
again positive-definite (see Add. |A"|) . 



Also for the present reduced description we can con- 
sider the semiclassical or Boltzmann limit; as described 
in Sect. Ill El this corresponds to neglecting the nondiag- 
onal matrix elements of the reduced density matrix, i.e., 



P{n ai }{n a2 } ~ f{n ai } 



S(„ o] }{n a2 } -Within such approxi- 
mation scheme the quantum-transport equation in i|59|) 
reduces to the following Boltzmann equation for the car- 
rier subsystem: 



df. 



{no} 



dt 



where 



Yl ( P {n a }{n a ,}f{n a ,} - P {n a ,}{n a }f{n a 



{»«'} 



(75) 



x ^ 2tt 

+ ~h 
q± 



N a 



1 1 
- ± - 

2 2 



-{n a } ~ £{n Q ,} 



±*0 



(76) 



r 



are the usual carrier-quasiparticle semiclassical scatter- 
ing rates given by the well-known Fermi's golden rule. 
We stress that, contrary to the global (carrier + quasi- 
particle) description considered in Sect. Ill El the scat- 
tering rates in (|76|l are not symmetric: P\ 



P 



{n a }{n a ,} 



this is a direct fingerprint of the irreversible 
nature of the transport problem induced by energy- 
relaxation/dephasing processes. 



G. Single-particle description 

Most of the electronic properties of interest in the anal- 
ysis of charge-transport phenomena in semiconductor 
nanostructures are single-particle quantities, i.e., physi- 
cal quantities ascribed to the generic particle in our elec- 
tronic subsystem, like carrier drift velocity mean kinetic 
energy, etc. In this case, the corresponding quantum- 
mechanical operator A is of the form: 



C a C a i 



(77) 



(«i = 0L2) correspond to the single-particle carrier distri- 
bution of the semiclassical Boltzmann theory, while the 
non diagonal contributions (pt\ ^ a-i) describe quantum- 
mechanical phase coherence between the single-particle 
states ai and a2- We stress that, while the reduced 
density-matrix operator p c describes the whole many- 
electron system, the single-particle operator p sp pro- 
vides an average or mean-field treatment of the carrier 
subsystem; indeed the latter fails in describing many- 
particle correlations, like Coulomb-correlation effects in 
quasi zero-dimensional systems^ 

Since p^ Q2 is the only quantity entering the evaluation 
of the average value in J7SJ, it is desirable to derive a 
corresponding equation of motion for the single-particle 
density- matrix in l|79l) : 



Jt' 



sp 



tr 



d 

Jt' 



(80) 



Inserting into the above expression the equation of mo- 
tion for the global density-matrix operator p in Ijl8|l we 
get: 



and its average value can be written as 



where 



Pa ia2 = tr {4 2 c ai p} 



(78) 



(79) 



is the so-called single-particle density matrix^ As we 
can see, this is defined as average of the product of cre- 
ation and destruction operators; its diagonal elements 



d 
Jt 

where 



— o* — — " 

dr aia2 dt' aia ' 



H o p 

h dt Paia ' 



-l — p s P 
c dt Paia '< 



dt 



H' 



, (81) 



(82) 



is the time variation induced by the noninteracting 
Hamiltonian H , 



dt 



o ' 

raic 



= —i tr 



{clc ai [n,P l (t)}} (83) 
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is the contribution due to the quantum-correlation oper- 
ator C in J22J), and 



dt 



p ' 



-tr {clc^iH, [It, p}]} 



(84) 



is the time evolution dictated by the scattering superop- 
crator T. 

For a better evaluation of the various contributions 
in (|82|) - (|84[1 it is convenient to expand the commutators 
entering the trace, regrouping the various terms in a dif- 
ferent way. More specifically, by inserting the explicit 
form of the single and double commutators, and using 
the cyclic property of the trace, we finally get 



— o ' 

dr aia2 



—p sp 

dV aia ' 



H, 



P 



-i tr • 



c t c 7t! 



'(*)} 



and 



— o ' 

dr aia ' 



-tr • 



H' 



c 1 c H 



X 



(85) 



(86) 



(87) 



As we can see, the various contributions to the time evo- 
lution of the single-particle density-matrix can be writ- 
ten as global average values of single as well as double- 
commutators. 

The term in (|85|l can be evaluated exactly. Indeed, re- 
calling the explicit form of the free-carrier + free-phonon 
Hamiltonian in @ and using the Fermionic anticommu- 

tation relations |c Q ,c^, | = S aa >, we obtain 



dt 



p sp 



C tt2 sp 



In contrast, for the first- and second-order interaction 
contributions in l(86|) and 1)87(1 it is not possible to ob- 
tain closed equations of motion for the single-particle 
density matrix p sp : indeed such contributions involve 
higher-order correlations, e.g., two-body and/or phonon- 
assisted density matrices^ In order to get a closed equa- 
tion for PcjZ a , an additional approximation is needed, 
the so-called mean-field approximation. The latter con- 
sists of a factorization of the higher-order correlation 
functions into products of single-particle density matri- 
ces p sp and/or quasiparticle populations Nq. The re- 
quired mean-field procedure and the explicit form of the 
resulting closed equation of motion depend on the partic- 
ular form of the interaction Hamiltonian considered, e.g., 
carrier-carrier, carrier-quasiparticle, etc. However, the 
free-evolution term l|85|) together with the first-order con- 
tribution in l|8()l) describe, in general, coherent phenom- 
ena — including Hartree-Fock renormalization and coher- 
ent phonons — while the second-order term in (|87|) de- 
scribe energy-relaxation/dephasing processes within the 
Markov approximation. 

At this point few comments are in order. The single- 
particle description discussed in this section is based on 



the Schrodinger picture: the equation of motion for p sp ia2 
[see Eq. H80fl ] is derived by treating the operators & and 
c in l|79l) as time-independent, while the time variation is 
fully attributed to the density-matrix operator p. Actu- 
ally, the most popular and commonly used approach 3 
to derive the equations of motion governing the time 
evolution of the single-particle density matrix is based 
on the Heisenberg picture: the density-matrix operator 
p entering Eq. I|79[l in the Heisenberg scheme is time- 
independent, while the time evolution is fully ascribed to 
the Fermionic operators via their corresponding Heisen- 
berg equations of motion: 



d - - 1 r* Tl 

TT" [Cot J £1 \ • 



dt 



(89) 



More precisely, within the Heisenberg picture Eq. I|8Q(I is 
replaced by: 



— s p - t / — 
dt Paia2 ~ r \dt 



itr{[ct 2 c QI ,F]p} (90) 



Contrary to the theoretical approach proposed in this 
paper, the usual Heisenberg treatment is based on a cor- 
relation expansion of the trace in l|9U|l : 3 starting again 
from the Hamiltonian separation in a hierarchy of 
kinetic equations involving higher-order density as well 
as quasiparticle-assisted density matrices is established; 
the different contributions are classified in terms of their 
perturbation order. Such infinite hierarchy is trun- 
cated/closed via the mean- field approximation previously 
recalled, and only at this level the Markov limit is usually 
introduced^ Aim of this paper, in contrast, is to analyze 
the Markov limit from a more general point of view; it is 
for this reason that the latter has been introduced in very 
general terms in Sect. Ill Cl before addressing any reduced 
description. 



1. Carrier-carrier interaction 

As first interaction mechanism we shall consider two- 
body Coulomb coupling. The corresponding interaction 
Hamiltonian can be written as 



V cc , , $ & c ,c , 

v cti cto iQi^a^ ^Qi ^2 a ± ot-2 



(91) 



where V cc , , is the Coulomb matrix element for the 

generic two-body transition a^a^ — > ol\Oli- 

In order to derive the explicit form of the second-order 
contribution to the single-particle dynamics in (|87() . two 
key quantities need to be evaluated: the inner commu- 
tator c]Lc Tit and the explicit form of the operator 

t. 

More specifically, by employing the anticommutation 
relations for the Fermionic operators, we get: 
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E 

030405 



V cc c 1 * c 1 c c 

0103,0405 02 03 04 05 



V cc & & c c ) 

' 0504,0302 ^a 5 ^o 4 ^03 "oi / 



(92) 



The above commutator has the same structure of the 
interaction Hamiltonian in (|91fl . i.e., it consists of a sum 
of products of four Fermionic operators. 



By inserting into Eq. Ij29(l the explicit form of the two- 
body Coulomb interaction Hamiltonian in (|91|) . we get: 



drU (r) I - 



E 



V cc , , $ $ c it 

y Q1Q2 ,0^2 a± &2 ck-l a 



(93) 



Recalling that 



^o(T)c^4 2 c a ,c a , Ul{T) = 44A;^e- 1 2 - Ql ° 2 ' , (94) 



we finally obtain 

^ / v 2 a i ft 2, ajaj ocia2,a' 1 a' 2 Qi «2 a x a 2 

(95) 

with 



OL\ Q.2 a 9 



1 

2^ 



t-to < e 



(96) 



We get again the same operatorial structure: a sum of 
products of four Fermionic operators. 

Given the two results in H92|) and H95|) , their commuta- 
tor — key ingredient in Eq. (|87|l — will involve, in general, 
products of six Fermionic operators. 

As anticipated, in order to get a closed equation of 
motion for the single-particle density matrix p sp we are 
forced to employ the mean-field approximation; the latter 
allows in this case to write, in general, the average values 
of six Fermionic operators as products of three single- 
particle density-matrix elements. 

By applying such mean-field factorization procedure 
to the explicit form of the outer commutator in l|87|l , the 
final result — not reported here — can be cast into the 
general form: 



dt 



raia-2 



r 



•cc,m r n sp~\ 



r 



cc, out \p S P 



(97) 



As for the case of the reduced description (see Sect. Ill ¥\ . 
the time variation of the single-particle density matrix 
is the sum of a positive — in-scattering — and a nega- 
tive — out-scattering — contribution. However, contrary 
to the global and reduced descriptions previously con- 
sidered, now the in- and out-scattering contributions are 
non-linear functions of the single-particle density matrix 
p sp . In particular, in this case of two-body interaction 
between a main (M) and a partner (P) carrier, the su- 
peroperators F cc > m and jF cc ' out both involve a product 



structure of the form p sp ™p spp (l - 



<)(X- 



Also for the present single-particle description it is pos- 
sible to consider the semiclassical or Boltzmann limit in- 
troduced in Sect. Ill El This amounts again to neglect 
non-diagonal density-matrix elements: p sp l0l2 = fS^S ai a 2 
. By inserting the above diagonal form of p sp into the in- 
and out-scattering superoperators f cc ' ln and J^ cc > aat } the 
following Boltzmann-like equation for the semiclassical 
single-particle distribution f^ p may be derived: 



dt 

where 



E K 1 - /o P )0/J - (1 - /SWa/oT 



(98) 
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o = f D 1 - K c wl 2 w« + ^- e «' - ^) (99) 

a a.' 



are two-body carrier-carrier scattering rates describing 
the main-carrier transition a! — > a accompanied by the 
partner-carrier transition a' — > a. As we can see, also in 
the semiclassical limit we deal with a non-linear transport 
equation; such nonlinearities are ascribed (i) to the pres- 
ence of the carrier distribution f sp of the initial partner 
carrier, and (ii) to the two Pauli-blocking factors (1 — f sp ) 
corresponding to the final states of both main and part- 
ner carriers. Comparing the semiclassical transport equa- 
tion in H98|) to its quantum-mechanical generalization in 
1|97JI. we clearly see that the various terms of the form 
(S aa i — p^ a , ) are the natural generalization of the Pauli- 
blocking factors (1 — f* p ) of the semiclassical theory. 



2. Carrier- quasiparticle interaction 

Let us now come to the carrier-quasiparticle coupling 
mechanism. By adopting as explicit form of the carrier- 
quasiparticle quantity 7i q = 7iL q the single-particle op- 
erator 



Tin 



h ^2 9aa ' 



q c a c a' ' 



(100) 



the carrier-quasiparticle interaction Hamiltonian in l|57f) 
is given by: 



aa' ,q 



where 



9aa',q = Sq/aa'.q = 5 q / 4>* a (r)e iqr Q , (r)rfr 



(101) 



(102) 



is the carrier-quasiparticle matrix element for the single- 
particle transition a' — > a induced by the quasiparticle 
bulk mode q. The explicit form of the coupling func- 
tion g depends on the particular carrier-quasiparticle in- 
teraction mechanism considered. In any case we have: 
5aa',q = 5q' q - q - As for carrier-carrier interaction, in or- 
der to derive the second-order contribution to the single- 
particle dynamics in (|87|l . we shall evaluate the inner 



commutator 



c H 



as well as the operator JC. 

More specifically, by employing again the anticommu- 
tation relations for the Fermionic operators we get: 



c f c K c 



-qp 



-y 



c 1 * tfc 



?Q 1 Q 3 ,qC Q , 2 & q C Q3 



(103) 



As we can see, the commutator has the same operatorial 
structure of the carrier-quasiparticle Hamiltonian, i.e., 
it involves products of one Bosonic and two Fermionic 
operators. 



By inserting into Eq. I|29[l the explicit form of the 
carrier-quasiparticle interaction Hamiltonian in (|lt)l(l , we 
get: 



9aa' ,qc]j^qC 



E/J(t) +h.c. 



(104) 



where h.c. denotes its hermitian conjugate. Recalling we finally obtain 
that 



2vr 



(105) 



jQC-qp _ _ 



eta.' ,q 



v c-qp, 



h.c. 



(106) 
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with 



where 



-qp,± 



1 

2^ 



t-t 



■dr . 



(107) 



We get again the same operatorial structure: one Bosonic 
times two Fermionic operators. 

Given the two results in ((103(1 and 1)106(1 , their commu- 
tator — key ingredient in Eq. 1)87(1 — will involve products 
of two Bosonic and two Fermionic operators as well as 
products of four Fermionic operators. As for the carrier- 
carrier interaction previously discussed, in order to get a 
closed equation of motion for the single-particle density 
matrix p sp we are forced to employ again the mean-field 
approximation; the latter allows to write (i) the aver- 
age value of two Fermionic times two Bosonic operators 
as the product of single-particle density-matrix elements 
Paa' times quasiparticle distributions, and (ii) the aver- 
age values of four Fermionic operators as products of two 
single-particle density-matrix elements. 

By applying such mean-field factorization procedure 
to the explicit form of the outer commutator in 1(87(1 , the 
final result not reported here — can be cast in the same 
form of the one for carrier-carrier interaction in ((97(1 . i.e., 



d_ 

dt f 



T 



-'Hi 



-<3P,outr spi 

[H \a\a 2 ■ 



(108) 

As for the case of carrier-carrier interaction, the above 
in- and out-scattering contributions are again non-linear 
functions of the single-particle density matrix p sp . More 
specifically, their general structure is of the form p sp (T — 
p sp ); Such nonlincarities — ascribed to Pauli-blocking 
effects — vanish in the so-called low-density limit (I — 
p sp 

In the semiclassical limit, by inserting the diagonal 
density-matrix form into the in- and out- scattering func- 
tional T c ~ qp ^ and ^~ip^\ W e ge t a non-linear Boltz- 
mann equation formally identical to the single-particle 
transport equation in 1(98(1 : 



dl 



■V 



dt 



E [(! - f« p )p c -j p fj a - oK> a 9p m , 

(109) 



pc-qp 

aa' 



T ^ l9ac 



' 8{e a 



±Cq) 



(110) 



q± 



denote semiclassical carrier-quasiparticle scattering rates 
for the single-particle transition a' — > a. We deal again 
with a non-linear transport equation; in this case, how- 
ever, such nonlincarities are only ascribed to the Pauli- 
blocking factor of the final state, and they vanish in the 
low-density limit as discussed in the following subsection. 



3. Low-density limit 



Let us finally consider the so-called "low-density limit" . 
To this aim, let us recall that within the single-particle 
description previously introduced the average occupation 
number for the generic state a is simply given by the 
diagonal elements of the single particle density matrix in 
®: 



/C = tr{at^} = /* 



(111) 



It is then clear that at low carrier concentrations (low 
densities) the magnitude of all density-matrix elements 



is much smaller than 1; More precisely, we have: p aa , 



sp 



0,(S o 



sp 
P aa' 



$ac 



In this limit the carrier-carrier scattering contribu- 
tions in 1(971) vanish since, as anticipated, they involve 
the single-particle density matrix of the partner carrier 
p spF . As a result, in the low-density limit the only 
non-vanishing contribution to the interaction dynamics 
is given by the carrier-quasiparticle terms in 1(108|) . More 
precisely, by inserting the low-density condition in our 
quantum transport equation 1(108(1 as well as in the cor- 
responding in- and out-scattering functions, we finally 
get the following linear transport equation: 



dV aia2 



sp 

a>2,u[a ! 2 P a ' a ' 



E 



r m , ,p 

a 1 a2,a' 1 a' 2 r, 



sp _ \ * pou 

a[a' 2 / y ai 



out 

a?.a\ 



sp 



(112) 
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We stress that the above linear in- and out-scattering op- 
erators are formally identical to the reduced-description 
ones in Eqs. (|73|l and l|74|l . This can be easily understood, 
considering that the present low-density limit is physi- 
cally equivalent to consider a system of just one electron 
interacting with the quasiparticle degrees of freedom; this 
would correspond to substitute the many-electron config- 
uration label {n a } with the state a of the only electron 
considered. 

Combining Eqs. (|88|l and (|112|) and neglecting the first- 
order term in (|86|) , we finally get the following quantum- 
transport equation: 



whose solution is simply given by: 



dt 



p sp 



(115) 



where 



^102 .a^ftj 



ih 



is the effective single-particle Liouville operator in the 
low-density limit. 



III. ANALYSIS OF PHASE-COHERENCE 
EFFECTS IN SEMICONDUCTOR 
NANOSTRUCTURES 

Let us now discuss in more detail the physical mean- 
ing of the in- and out-scattering operators in Eqs. (|113|) - 
(|114ll . As shown in Eq. I|112|) . the generic matrix ele- 
ment r in / out describes how the density-matrix ele- 

a±a 2 ,a 1 a 2 J 

ment p a > a i influences the time evolution of the clement 
p ai a 2 - We shall try to classify such matrix elements ac- 
cording to their role and physical meaning. 

As a first class of terms, let us consider the semiclassi- 
cal ones (aia^ = a 2 a' 2 ): 



Ti 



r 



act. a' a' 



Pr, 



/ - 1 a" a' 



(117) 



As discussed in App. [XJ the latter, also referred to as T\ 
contributions, are fully expressed in terms of the semi- 
classical rates P aa > [see Eq. (jllUf) ] and act on the carrier 
distribution f a = p a p a only giving rise to electron and 
energy redistribution within the traditional Boltzmann 
picture. 

A second class of terms is given by the so-called T 2 
contributions (a\a 2 = a^a^): 

la\a 2 = ^ a 1 a 2 ,a 1 a 2 = — I y ^ P a " ai + ^ ] Pg"a 2 J ■ 

\ a" a" ) 

(118) 

In the absence of other T matrix elements, these T 2 con- 
tributions describe the damping of non-diagonal density- 
matrix elements p ai a 2 — PZ^a^ a ^ so cane d inter-state 
polarizations; 3 Indeed, from Eq. (|1 L2|) we easily get: 

d 



dt 



Pa±c 



c-qp 



T 2 



(119) 



v T 2 



p aia2 (t) = e^2^o) paia2 ^ . ( 120 ) 

The above exponential decay (7^ a2 1S a l wa Y s negative) 
of non-diagonal density-matrix elements or inter-state 
polarizations is the well-known decoherence or dephas- 
ing process £ We stress that, in spite of the quantum- 
mechanical (i.e. coherent) character of p ai a 2 , the de- 
phasing rates j T2 in 1118fl involve semiclassical ingredi- 
ents only; indeed, the latter can be regarded as the sum 
of the two semiclassical out-scattering rates for states a\ 
and a 2 (see Add. lA^I . 

In addition to these two classes of terms, acting inde- 
pendently on f a and p ai a 2 , the fully non-diagonal de- 
scription in 1112fl involves also non trivial coupling con- 
tributions, namely (i) coupling between different polar- 
izations (p aia2 — * Pa[a' 2 )i an d (h) terms of the form 
aa — > a±a 2 and vice versa; the latter describe coupling 
between f a and p ai a 2 - Denoting with T^-contributions 
these last / — > p terms and neglecting p — > p' couplings, 
the quantum-transport equation in (|112f> can be schemat- 
ically written as: 



d_ 

dt 



7 Tl 7 



T 3 



(121) 
p coupling 



where j T3 denotes schematically all the / 
terms previously mentioned. 

We clearly see that, in the absence of X3 terms (7 T3 = 
0) there is absolutely no coupling between population 
(/) and polarization (p). This approximation scheme, 
known as T1T2 model* 3 - is the most popular model used 
for the interpretation of ultrafast optical experiments in 
semiconductors . 

In contrast, the introduction of these T3 terms — not 
included in the conventional T1T2 model — gives rise to 
a non-trivial coupling between / and p, which manifests 
itself in a residual single-particle phase coherence also in 
steady-state conditions. More specifically, while within 
the T1T2 model in the limit t — ► 00 all polarization terms 
p ai a 2 vanish [see Eq. (|120fl ] and the resulting steady- 
state density matrix is diagonal {p sp ia2 — f ai S ai a 2 ), in 
the presence of T3 terms we may have a steady-state so- 
lution characterized by non-zero values of p ai a 2 ( see be- 
low). As discussed in App.^] in the small-coupling limit 
the steady-state solution p a p ±a2 of the transport equation 
in 11151) is always positive-definite. This suggests the in- 
troduction of the unitary transformation T aa = {q\q) 
which diagonalizes the steady-state solution p sp . This 
allows to clearly express the positive character of the so- 
lution according to Eq. (|38f) . 

In order to gain more insight into the fully non- 
diagonal density-matrix approach presented so far, let us 
consider an extremely simplified scenario: an electronic 
two-level system interacting with a single phonon mode 
q. In this case we deal with a two-by-two density matrix 
of the form: 



(122) 



pZ 


pZ) 






[pZ 


pZ) 


={ f ; 


h) 
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Here, the diagonal elements p a p a and p s bb coincide with 
the semiclassical ground- and excited-state level popula- 
tions f a and ft,, while the non diagonal element p = p s b p a 
(together with its complex conjugate p* — p sp b ) describes 
the degree of quantum-mechanical phase coherence be- 
tween states/levels a and b. Let us introduce the follow- 
ing (arbitrary) v — {a\,a2\ mapping: 1 = {a, a}, 2 = 
{6, 6}, 3 = {b, a}, 4 = {a, b}. Within such representa- 
tion, the two-by-two density matrix in l|122|) is mapped 
into a four-dimensional vector, and the single-particle 
Liouville superoperator £ in (|116[) will correspond to a 



four-by-four matrix. More specifically, within the four- 
dimensional mapping given above the transport equation 
pi5|l in steady-state conditions reduces to the following 
homogeneous linear problem: 



£aa,aa 


£aa,bb 


£aa,ba 


£aa,ab 






£bb,aa 


£bb,bb 


£bb,ba 


£>bb,ab 






£ba,aa 


£ba,bb 


£-ba,ba 


£-ba,ab 






£ab,aa 


£ab,bb 


£-ab,ba 


£-ab,ab J 





= (123) 



with 
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(124) 



Here Ae = e& — e a denotes the interlevel energy splitting. 
Starting from the explicit form of the electron-quasipar- 
ticle in- and out-scattering superoperators in I|113I114|) 
and assuming real and symmetric single-particle wave- 



functions 4>ot (r) in Eq. (|102|l , the resulting expression for 
the matrix elements of T for our two-level system is of 
the form: 



^bb,aa 



aa,bb 


Taa,6a 


^aa,ab 


\b,bb 


Fbb,ba 


^bb,ab 


ba,bb 


±ba,ba 


^ba,ab 


ab,bb 


r ^ ab.ba 


r ^ ab.ab 




(N + l) 
-(N + l) 

K(^v + i) 



o 



-i(2JV + l) 
|(2JV + 1) 




(125) 



where £ = ^''-^^/"^ (with e<j = Ae), e c P denotes 
the carrier-phonon coupling constant, and N is the Bose 
occupation number in l|61|l for our two-level system. 

As already stressed, by neglecting / — > p terms (which 
corresponds to set g aa ,q = fff>b,q)> populations and polar- 
izations are totally decoupled. In contrast, in the pres- 
ence of these T 3 terms our steady-state solution exhibits 
a residual single-particle coherence, i.e., a polarization p 
different from zero. 

Figure ^ shows the modulus of the interlevel polariza- 
tion p as a function of the dimensionless coupling pa- 
rameter a = In this numerical example we have 
chosen Ae = 25meV, T = 300 K, and a = b = 0.25. 
As we can see, \p\ comes out to be proportional to the 
coupling parameter er, as can be readily verified by a 
closer inspection of our four-by-four superoperator C in 
l)124[l . It follows that for any finite value of the coupling 
parameter a, we deal with a non-diagonal steady-state 
solution p [see Eq. (|122f) ]. More specifically, its diagonal 
elements f a and /& are cr-independent and correspond 
to the semiclassical equilibrium distribution: fy~ = jV ^ 1 ; 



the non-diagonal elements p and p* will increase linearly 
with the coupling parameter a. In the small-coupling 
limit (cr < 1) such non-diagonal contributions are much 
smaller than the diagonal ones. For increasing values of a 
we progressively enter the strong-coupling regime, up to 
the point where our two-by-two density matrix in l|122|) 
is no more positive definite. Indeed, for \p\ > \fjah, 
the determinant of p (i.e., f a fb — \p\ 2 ) becomes negative. 
As anticipated, this is exactly the regime for which the 
Markov approximation itself is no longer valid. How- 
ever, for small and moderate values of a our two-by-two 
density matrix is positive-definite, which suggests the in- 
troduction of a dressed-state basis in which the latter is 
diagonal. The new populations f a and f b (dashed curves 
in Fig. ^) can be regarded as the average occupation of 
such dressed states. As we can see, for a = they co- 
incide with the noninteracting thermal ones; for increas- 
ing values of the carrier-phonon coupling the population 

ratio =£■ decreases. Such a behavior can be physically 

described in terms of a phonon-induced renormalization 
of the interlevel energy splitting Ae. Indeed, such renor- 
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FIG. 1: Phonon-induced single-particle phase coherence foi 
an isolated quantum-dot system: modulus of the interleve 
polarization \p\ (solid curve) and dressed-states populations /; 
and f a (dashed curves) as a function of the coupling-constanl 
ratio cr, for an interlevel energy splitting Ae = 25 meV a1 
room temperature. 



malized transition energy can be also obtained from the 
imaginary parts of the non-zero eigenvalues of the effec- 
tive Liouville operator in l|124|) . 

We finally stress that such bidimensional dressed basif 
as well as the corresponding energy shift can be regardec 
as the simplest example of polaronic phase-coherence anc 
state renormalization (see below). 



IV. POLARONIC COHERENCE IN BIASED 
SEMICONDUCTOR SUPERLATTICES 

As a concrete example of scattering-induced single- 
particle phase coherence in state-of-the-art semiconduc- 
tor nanostructures, let us consider high-field transport in 
biased semiconductor superlattices. As extensively dis- 
cussed in Ref. l27l the treatment of carrier-phonon scat- 
tering in the presence of strong electric fields requires a 
gauge-invariant formulation of the problem; the density- 
matrix treatment proposed in Ref. [27] has shown that the 
only single-particle basis a in which the Markov limit is 
properly defined is the Wannier-Stark one. 

In order to gain more insight into the single-particle 
phase coherence previously discussed, we have performed 
fully three-dimensional calculations of high-field charge 
transport in state-of-the-art semiconductor nanostruc- 
tures. In particular, aim of our description was to prop- 
erly treat — in addition to the carrier quantum confine- 
ment along the growth direction — the in-plane energy- 
relaxation and thermalization dynamics. 

As prototypical system, we consider a state-of-the- 
art GaAs-based nanometric superlattice. More specifi- 
cally, we have performed a detailed investigation of the 
30A GaAs/Al .3Gao.7As superlattice structure shown in 
Fig. |21 The single-particle carrier states {\a}} are de- 
scribed within the usual envelope-function approxima- 
tion in terms of a space-independent effective mass m*. 




z(nm) 



FIG. 2: Schematics of the prototypical 30A/30A 

GaAs/Alo.3Gao.7As superlattice structure considered in our 
simulated experiments: Real-space periodic nanostructure 
profile (shaded regions correspond to barriers) and charge dis- 
tribution corresponding to the ground-state (k« = 0) envelope 
function in the field-free case[(n(r||) oc \<$, =0iI/=1 (r||)| 2 )]. 

They come out to be products of two-dimensional plane 
waves and one-dimensional envelope functions: 

^(r) = -^e lk - r ^ll || (r||) , (126) 

n denoting a suitable normalization area. 

In the free-field case, the envelope functions 0" in l|126|) 
reduce to one-dimensional Bloch states 0| v correspond- 
ing to the periodic heterostructure potential reported in 
the inset of Fig. [21 As we can see, we deal with a rel- 
atively small band-edge discontinuity (V = 250 meV). 
The latter, combined with a barrier width of 30 A gives 
rise to significant inter- well carrier tunnelling. This is 
confirmed by the field-free ground-state charge distribu- 
tion (solid curve in Fig. [2J, which shows a clear finger- 
print of carrier derealization. 

The inter-well coupling displayed in Fig. [21 should 
translate into a dispersive energy-momentum relation 
along the growth direction. This is confirmed by the 
superlattice miniband profile reported in Fig. [31 As 
we can see, we deal with a carrier miniband only: Its 
width amounts to about 90 meV and it is therefore larger 
than the longitudinal optical (LO)-phonon energy (~ 36 
meV). 

In the presence of an applied field F along the growth 
direction, the one-dimensional envelope functions 0" 
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FIG. 3: Single-miniband diagram (miniband width of about 
20 meV) corresponding to the superlattice structure depicted 
in Fig. |21 The superlattice potential profile (band-edge dis- 
continuity of 0.25 eV) is also shown in the inset. 

within the scalar-potential gauge correspond to the usual 
Wannier-Stark states^ The latter are displayed in Fig.0] 
for different values of F. As we can see, for increasing 
values of the field we deal with an increasing state lo- 
calization and a corresponding suppression of inter-well 
single-particle tunnelling. 

Primary goal of our simulated experiments was the 
study of the current-voltage characteristics of the super- 
lattice structure previously introduced, in the presence 
of carrier-LO phonon scattering. In particular, a thermal 
bath of three-dimensional (bulk) dispersionless longitudi- 
nal polar optical phonons has been assumed, employing 
the carrier-phonon interaction model given in Ref. l28t 
We have therefore evaluated the carrier drift velocity as 
a function of the applied field. The latter can be readily 
computed according to Eq. l(75|) . using as single particle 
quantity the velocity operator: 

P ih 

A sp = v = = V . (127) 

TO* TO* 

More specifically, we get: 

« d = £w^ Q > (128) 

aoc f 

where p sp a , is the single-particle density matrix, and v aa > 



FIG. 4: Charge distribution corresponding to the Wannier- 
Stark states in the superlattice structure of Fig.|5|for different 
values of the applied field F. For each field three states are 
displayed: n = — 1 (dotted curve), n = (solid curve), and 
n = 1 (dash-dotted curve). 

denote the matrix elements of the velocity operator l|127|) 
within our a representation. 

FigureElshows the steady-state carrier drift velocity as 
a function of the applied field for the superlattice struc- 
ture of Fig. [21 at room temperature in the low-density 
limit. The peak at ~ 70 kV/cm corresponds to the 
phonon resonance, i.e., for this value of the applied field, 
the Wannier-Stark energy eFd — d being the superlattice 
period — is equal to the LO-phonon energy. We stress 
that such feature well agrees with the results obtained 
with the non-equilibrium Green's function treatment by 
Jauho and co-workers ^ As we can see, we get drift- 
velocity values significantly different from zero. This re- 
mark is particularly important in view of the fact that 
such field-induced current is entirely due to the pres- 
ence of non-diagonal density-matrix elements; indeed, 
within the semiclassical limit previously introduced we 
have p sp a , = fa P 8aa>, and therefore Eq. (|128|l reduces to: 

V d = J2 V <*«fa ■ ( 129 ) 

a 
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FIG. 5: Steady-state carrier drift velocity as a function of the 
applied field for the superlattice structure of Fig. at room 
temperature, in the low-density limit. The line is a guide to 
the eye. 



As we can see, the only matrix elements of the velocity 
operator involved are the diagonal ones, i.e., v a=a >\ for 
localized as well as periodic basis states a such diago- 
nal matrix elements are known to vanish, and therefore 
within the semiclassical limit the drift velocity in l|128|) 
is equal to zero, i.e., no current. 

It is then easy to conclude that the significantly large 
values of the drift velocity reported in Fig. are en- 
tirely due to the non-diagonal density-matrix elements 
Pa^a' entering Eq. (|128f) : in turn such non-diagonal con- 
tributions reflect a scattering induced phase coherence 
between our Wannier-Stark states. As a result, we can 
say that the significant carrier drift velocity in Fig. [S] 
is an unambiguous fingerprint of phonon-induced phase 
coherence in semiconductor superlattices. 



that the Markov limit alone — i.e., not combined with the 
diagonal or semiclassical approximation — may lead to 
totally non-physical results: (i) within the global (carrier 
+ quasiparticle) picture the Markov limit introduces a 
fictitious energy redistribution/relaxation and dephasing 
which has no physical counterpart within a closed-system 
scenario; (ii) both for the global and for the reduced pic- 
tures, the resulting effective Liouville superoperators do 
not correspond, in general, to a CP map, thus prevent- 
ing from any "robust" time-dependent solution of our 
quantum-transport equation. 

A second important remark is that the combination 
of the Markov limit with the semiclassical approxima- 
tion — i.e., the neglect of non-diagonal density-matrix 
elements — leads, also for the global description, to 
Boltzmann-like equations whose solution can be shown 
to be positive-definite at any time. 

In spite of the fact that our density-matrix formu- 
lation within the usual Markov limit does not trans- 
late into CP maps, our analysis has shown that within 
the small-coupling limit the proposed quantum-transport 
equations are always characterized by a positive-definite 
steady-state solution. 
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APPENDIX A: GENERAL PROPERTIES OF THE 
SCATTERING SUPEROPERATOR 



V. SUMMARY AND CONCLUSIONS 

We have presented a general density-matrix formula- 
tion of quantum-transport phenomena in semiconductor 
nanostructures typically employed in the design of state- 
of-the-art quantum devices. More specifically, contrary 
to the conventional single-particle correlation expansion, 
we have investigated separately the effects of the Markov 
limit and of the reduction procedure. Our fully operato- 
rial approach has allowed us to better identify the general 
properties of the scattering superoperators entering our 
effective quantum-transport theory at various description 
levels, e.g., N electrons-plus-quasiparticles, TV electrons 
only, and single-particle picture. 

In addition to coherent transport phenomena charac- 
terizing the transient response of the system, the pro- 
posed theoretical description has allowed the study of 
scattering induced phase coherence in steady-state con- 
ditions. In particular, based on the proposed approach 
we have investigated polaronic effects in strongly biased 
semiconductor superlattices. 

One important conclusion of our theoretical analysis is 



Aim of this appendix is to discuss the general structure 
and symmetry properties which characterize the global 
(carrier + quasiparticle) in- and out-scattering superop- 
erators in (|33() - (|34|l . their reduced versions in Q73 [l -I|74 [l . 
and the single-particle ones in 1)1 13(1 - (|114JI . 

To this end, let us consider first the semiclassical limit 
introduced in Sect. Ill El For the global picture as well as 
for the reduced and single-particle descriptions we deal 
with a Boltzmann transport equation of the form: 

i' 

In analogy with the definition of the scattering super- 
operator T in lt2*rj|) . the transport equation (|A1|) can be 
rewritten as: 

§=Z>*'' ( A2 ) 

i' 

where 

7»' = 1& - 7^ = Pw - S W J2 P"* ( A3 ) 

i" 
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is again the difference of in- and out-contributions. As a 
first remark, we notice that (i) the in-contribution coin- 
cides with the semiclassical scattering rate Pu> and (ii) 
the out-contribution is a function of the in-one: 

i" 

The Boltzmann transport equation in (|A1|) is supposed 
to preserve the total number of particles N — /,;, i.e., 

f = Ef = I>/, = o. w 

i ii' 

The above conservation law — which corresponds to 
the trace-preserving property of the quantum-mechanical 
case — needs to be verified for any carrier distribution /,•; 
this requires that: 

£>«'=0, (A6) 

i 

which can be readily verified starting from the definition 
of the operator 7 in ljA3fl . 

Let us now consider the eigenvalue problem for the 
semiclassical scattering operator 7: 

X>'/i'=A/<. (A7) 

i' 

As a first property, the eigenvalue spectrum is non- 
positive definite: A < 0; secondly, the particle-preserving 
property in (|A6|) requires the existence of the A = 
eigenvalue, whose eigenvector fa corresponds to the 
steady-state solution of the problem. 



We shall finally discuss the positive-definite character 
of the semiclassical carrier dynamics. Starting from the 
Boltzmann equation in l|A2|l . we can express the carrier 
distribution at time t in terms of the semiclassical prop- 
agator Sij (t — to). 

/<(t)=£*e(i-*o)/i(to), (A8) 
3 

where Si?(0) = Sij and 

^ = Ewi' (A9) 

i' 

By inserting into the above equation the explicit form of 
the scattering operator 7 in l|A3J) we get: 

=-7i*tf , (A10) 
v 

where % = P,/j denotes the total or out-scattering 
rate for state i. The above differential equation can be 
formally integrated according to: 

s ij (t)=e-^ t - t ^S ij + [ dt'e-^-^YsPii'Si'M) ■ 

J to v 

(All) 

This formal solution — also known as Chamber's for- 
mulation^ — is the starting point of the semiclassical 
Neumann expansion. Indeed by iteratively substituting 
Eq. (| All|> into itself, we obtain: 



;(*) = £/ dti.../ dt n e-^(^)p iil e-7n(* 1 -* 2 )... Pi _ i . e -7,(*n-ta) , 

— n tn *° u,...,i„_i 



(A12) 



n=0 



The above expansion can be regarded as the semiclassical 
counterpart of the Neumann series in Ijl5[) i2fli By recalling 
that the scattering rates Pw and j i are positive-definite 
quantities, it is easy to recognize that all the terms of the 
above Neumann expansion are non-negative (i.e., greater 
or equal to zero); therefore, we can conclude that all the 
matrix elements s^- (t) of the semiclassical propagator are 
positive-definite, which in turn ensures that the distri- 
bution function /j in i|A8(l is itself positive-definite. As 
anticipated, we see that the combination of the Markov 
limit with the semiclassical or diagonal approximation 
leads to a completely positive time evolution, both for 
the global (carrier + quasiparticle) case and for the re- 
duced as well as single-particle pictures. 



Let us now move to the quantum-mechanical case. A 
closer inspection of the in- and out-scattering superoper- 
ators in (|33|> - (|34|> as well as of their reduced versions in 
!|73 )l -l|71 )l and (|113I) - H114|) reveals that they have a general 
structure of the form: 

r in . ., ., = V in - ■, ■, + V m * •/ •/ (A13) 

*l*25*i* 2 *1*25* 1 * 2 *2*1 5 * 2 Z 1 ' \ ^ J 

pout _ pout 1 pout* (A 141 

*l*25*i* 2 *1*2 ) * 1 * 2 *2*1;* 2 Z 1 ' J 

with 

Ki 2 ,i'A = ^9ni>9ii'V* 2i ' 2 (A15) 
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and 



-pout 

i\i 2 ,i'-yi' 2 



*2 l 2 / j ,%\l' x 



(A16) 



Exactly as in the scmiclassical theory previously recalled 
[see Eq. (|A4(l j. the out-scattering superoperator is fully 
determined by the in-one. The latter, in turn, can be 
regarded as the quantum-mechanical (i.e., nondiagonal) 
generalization of the conventional Fermi's golden rule; 
indeed, its semiclassical contributions = £2^2) are °f 
the form: 



pin 



2vr 



\gu>\ (Pa'+Vi'i) 



(A17) 



Similar to the particle-conservation law in 1A5(I . in the 
quantum-mechanical case we ask the scattering superop- 
erator 



p -pin pout 

Hi2,i'-ii' 2 i\i 2 S' x i' 2 hi2,i' 1 i' 2 



to be trace preserving, i.e., 



-tr{p}=tr{r(p)} = 



(A18) 



(A19) 



Rewriting the above equation within our generic i-basis 
we get: 



EE r 







(A20) 



This relation should hold for any generic p; this requires 
that: 







(A21) 



The above trace-preserving condition can be readily veri- 
fied starting from the general structure of the in- and out- 
scattering superoperators in (|A13|) - (|A14|) . We stress that 
such property is verified for any form of the in-scattering 
superoperator P 1 ", i.e., does not depend on the specific 
structure of the scattering operator in (|A15|) . 

Let us now come to the eigenvalue problem for the scat- 
tering superoperator T and for the corresponding effec- 
tive Liouville superoperator L [see Eqs. (|115ll and l|116l) ]. 
By denoting with j = i\i2 the generic pair of indices for 
our density matrix, we have: 



(A22) 



Contrary to the semiclassical picture, now the eigenval- 
ues A are complex quantities and, more important their 
real parts may assume positive as well as negative val- 
ues. This is a clear indication of the fact that we are 
not dealing with a CP map (see below). However, the 
trace-preserving property in l|A21(l requires also for the 
quantum-mechanical case the existence of the A = 
eigenvalue, whose eigenvector pj = pi x i 2 corresponds to 
the steady-state solution of the problem. 



We shall finally discuss the possible positive-definite 
character of our density-matrix p. In very general terms 
we can express the density matrix at time t via a 
quantum-mechanical propagator or quantum-mechanical 
map S: 

Phi 2 (t) = y^^i»2,jij 2 (* - t o)Pj 1 j 2 {to) , (A23) 
where 5 JlJ2jlj2 (0) = S ili2tjlj2 and 



^_Siii 2 ,jij 2 — ^ ( £i\i2,i' 1 i' 2 Si' x i 2 ,ixj 2 



(A24) 



Using the more compact operatorial notation we have: 
p(t) = S (p(t )) (A25) 

with 

dS 



dt 



cs . 



(A26) 



It is well known that the most general form of a CP map 
is given by: 

S(p(t )) = Y / Mk(t-t )p(t )MHt-t Q ) , (A27) 



where the generic set of operators A4k — usually referred 
to as "Krauss operators" — should obey the normaliza- 
tion condition 



E 



MkMl 



J . 



(A28) 



A particular class of CP maps is given by the so-called 
Lindblad maps, whose effective Liouville superoperators 
are of the form: 



C{p) = hH a ,p] 
in 



E 



A, 



.Ml 



(A29) 

For the particular case of Hermitian operators Ak, the 
Lindblad form in Eq. (|24|l is recovered. 

Let us now try to compare the effective scattering dy- 
namics described by the transport equation in l|115|) with 
the general Krauss expansion in l|A27|) and with the Lind- 
blad form in JX29J. 

For the global description, the double-commutator in 
(|19fl is not Lindblad-like; indeed, as already pointed 
out, its eigenvalue spectrum is in general not negative- 
definite. This tells us that, while the steady-state solu- 
tion of the Lindblad form in H24fl is again proportional to 
the identity operator I, its dynamical evolution in gen- 
eral does not preserve the positive-definite character of p 
and may also exhibit singularities. Such anomalous and 
totally non-physical behavior — not present within the 
semiclassical or Boltzmann picture — is again the result 
of the Markov limit. 
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Also for the case of the reduced and single-particle de- 
scriptions, their effective Liouville operators derived in 
Sects. Ill Fl and Iff Gl do not correspond, in general, to 
CP maps, and its time evolution does not preserve the 
positive-definite character of p. 

In spite of the fact that, in general, we are not deal- 
ing with CP maps, it is possible to show that for all 



the transport equations derived in the paper we have a 
steady-state solution. Moreover, in the so-called small- 
coupling limit, i.e., when the perturbation Hamiltonian 
H' is small compared to the nonintcracting Hamiltonian 
H , it is possible to show that such steady-state (A = 0) 
solution is always positive-definite [see Eq. (|38Jl ] . 
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